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1. Introduction 

The jet correlations are interesting in the context of recent detailed studies of hadron-hadron 
correlations in nucleus-nucleus ^ and proton-proton collisions. Those studies provide inter- 
esting information on the dynamics of nuclear and elementary collisions. Effects of geometrical 
jet structure were discussed recently in Ref.[^. No QCD calculation of parton radiation was per- 
formed up to now in this context. Before going into hadron-hadron correlations it seems indis- 
pensable to better understand correlations between jets due to the QCD radiation. In this paper 
we address the case of elementary hadronic collisions in order to avoid complicated and not yet 
well understood nuclear effects. Our analysis should be considered as a first step in order to un- 
derstand the nuclear case in the future. In leading-order coUinear-factorization approach jets are 
produced back-to-back. These leading-order jets are therefore not included into correlation func- 
tion, although they contribute a big (~ ^) fraction to the inclusive cross section. 

The truly internal momentum distribution of partons in hadrons due to Fermi motion (usually 
neglected in the literature) and/or any soft emission would lead to a deconelation from the simple 
kinematical back-to-back configuration. In the fixed-order coUinear approach only next-to-leading 
order terms lead to nonvanishing cross sections at / ;r and/or / P2,t (moduh of transverse 
momenta of outgoing partons). In the /Cf -factorization approach, where transverse momenta of 
gluons entering the hard process are included explicitly, the decorrelations come naturally in a 
relatively easy to calculate way. In Fig.|I] we show /c, -factorization processes discussed up to now 
in the literature [Q, |^ ^]. The soft emissions, not explicit in our calculation, are hidden in model 
unintegrated gluon distribution functions (UGDF). In our calculation the last objects are assumed 
to be given and are taken from the literature [||, |9|,p^. 



emissions 




Figure 1: Diagrams for -factorization approach included in the literature. We shall call them A\ and A2 
for brevity. 

In addition we include two new processes (see Fig.|2|), not discussed up to now in the context 
of ^, -factorization approach. We shall discuss their role at the RHIC energy W = 200 GeV. 

Furthermore we compare results obtained within the -factorization approach and results ob- 
tained in the NLO coUinear-factorization. Here we wish to address the problem of the relation 
between both approaches. We shall identify the regions of the phase space where the hard 2^3 



2 



Dijet correlations in pp collisions at RHIC 



Antoni Szczurek 




Figure 2: Diagrams for ^, -factorization approach included for the first time here. We shall call them Bi and 
B2 for brevity. 

processes, not explicitly included in the leading-order A;f -factorization approach, dominate over the 
2 — > 2 contributions calculated with UGDFs. We shall show how the results depend on UGDFs 
used. 

We shall concentrate on the region of relatively semi-hard jets, i.e. on the region related to the 
recently measured hadron-hadron correlations at RHIC. 

2. Formalism 

The cross section for the production of a pair of partons (k,l) can be written as 

da{hih2^ jetjet) f , , cfku (f-k2j 1 , 77772 

5^{tl,t + 1^2,,-tl,t-y2.t)-^i{^lA.t)^A^2A,,) , (2-1) 



where 



^^ = !!!Lie+^'+^^e+^^ (2.2) 
X2 = • ' + ' ' , (2.3) 



and m\^t and m2,t are so-called transverse masses defined as m,i = y pj^ + rrr, where m is the mass 
of a parton. In the following we shall assume that all partons are massless. The objects denoted 



by ^i{xi^k\j) and ^j{x2-,k2,) in the equation (gj) above are the unintegrated parton distributions 



in hadron h\ and /j2, respectively. They are functions of longitudinal momentum fraction and 
transverse momentum of the incoming (virtual) parton. 

After some simple algebra one obtains a handy formula: 

da{pi,t,P2.t) 1 1 



dp\fdp2,t 2 



• ^ - 471 / p\,iP2.t XI / '^y^'^y'i- \qtdqtd(^q, 

i,j,U 

( 16k^{x,X2sY ^-^^'^ " kl)\^^tix,,kl)^j{x2,kl)^ . (2.4) 
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This 5 -dimensional integral is now calculated for each point on the map pi t x p2^f 

Up to now we have considered only processes with two explicit hard partons. Now we shall 
discuss also processes with three explicit hard partons. In Fig.^ we show a typical 2 —>■ 3 process 
and kinematical variables needed in the description of the process. We select the particle 1 and 
2 as those which correlations are studied. This is only formal as all possible combinations are 
considered in real calculations. 




Figure 3: A typical 2^3 process. The kinematical variables used are shown explicitly. 

The cross section for /ji/j2 — > gggX can be calculated according to the standard parton model 
formula: 

da{hih2^ ggg) = J dxidx2g\{x\,ix^)g2{x2,ix^) dd{gg ^ ggg) , (2.5) 
where the longitudinal momentum fractions are evaluated as 

xi = ^exp(+3;i) + — ^rexp(+3;2) + ^exp(+3;3) , 

^yS ^/S ^/ S 

X2 = -^exp{-yi) + ^exp{-y2) + ^exp{-y3) . (2.6) 
^y S S y/ s 

After a simple algebra |Q] we get finally: 

1 



da = 54^4^2 ^igi(^bM/)-^2g2(-^2,M/) \^2^?,\'^P\.tdp\.tP2,tdp2,td<^-dy\dy2dy^ , (2.7) 

where 0_ is restricted to the interval (0, n). The last formula is very useful in calculating the cross 
section for particle 1 and particle 2 correlations. 

3. Results 

Let us concentrate first on 2 ^ 2 processes calculated within fc, -factorization approach. 

In Fig.0 we show two-dimensional maps in {pij,p2.t) for all ^, -factorization processes shown 
m Fig.| and Fig.|[ Only very few approaches in the literature include both gluons and quarks and 
antiquarks. In the calculation above we have used Kwieciriski UPDFs with exponential nonper- 
turbative form factor (bo = 1 GeV^') and the factorization scale = [pt^min + Pt,maxY = 100 
GeV2. 

In Fig]5| we show fractional contributions (individual component to the sum of all four compo- 
nents) of the above four processes on the two-dimensional map (3'i,3'2)- One point here requires a 
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Figure 4: Two-dimensional distributions in pi j and p2.t for different subprocesses gg gg (left upper) 
88 ~^ 1^ (right upper), gq gq (left lower) and qg qg (right lower). In this calculation W = 200 GeV and 
Kwiecihski UPDFs with exponential nonperturbative form factor (/jq = 1 GeV^^) and /i^ = 100 GeV^ were 
used. Here integration over full range of parton rapidities was made. 

better clarification. Experimentally it is not possible to distinguish gluon and quark/antiquark jets. 
Therefore in our calculation of the {y\,y2) dependence one has to symmetrize the cross section 
(not the amplitude) with respect to gluon - quark/antiquark exchange {y\ — > 3^2, 3^2 Ji)- While 
at midrapidities the contribution of diagram B\ + 62 is comparable to the diagram Ai, at larger 
rapidities the contributions of diagrams of the type B dominate. The contribution of diagram A2 is 
relatively small in the whole phase space. When calculating the contributions of the diagram A\ 
and A2 one has to be careful about coUinear singularity which leads to a significant enhancement 
of the cross section at 0_=O and y\ = j2, i.e. in the one jet case, when both partons are emitted 
in the same direction. This is particularly important for the matrix elements obtained by the naive 
analytic continuation from the formula for on-shell initial partons. The effect can be, however, 
easily eliminated with the jet-cone separation algorithm 

For completeness in Fig.^ we show azimuthal angle dependence of the cross section for all 
four components. There is no sizeable difference in the shape of azimuthal distribution for different 
components. 

The Kwieciilski approach allows to separate the unknown perturbative effects incorporated via 
nonperturbative form factors and the genuine effects of QCD evolution. The Kwiecinski distribu- 
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Figure 5: Two-dimensional distributions of fractional contributions of different subprocesses as a function 
of yi and y2 for gg gg (left upper) gg — > qq (right upper), gq gq (left lower) and qg — > qg (right lower). 
In this calculation W = 200 GeV and Kwieciriski UPDFs with exponential nonperturbative form factor and 
bo = 1 GeV^' were used. The integration is made for jets from the transverse momentum interval: 5 GeV 

<pi^t,P2,t <20 GeV. 

tions have two external parameters: 

• the parameter bo responsible for nonperturbative effect (for details see [Q]), 

• the evolution scale /i^ (for details see [Q]). 

While the latter can be identified physically with characteristic kinematical quantities in the pro- 
cess jU^ ~ plt^Pit^ ^^'^ ^^^^ ^'^^ °f nonperturbative origin and cannot be calculated from first 
principles. The shapes of distributions depends, however, strongly on the value of the parameter 
bQ in which the inital momentum distribution is encoded. This is demonstrated in Fig.^ where 
we show angular correlations in azimuth for the gg — > gg subprocess. The smaller bo the bigger 
decorrelation in azimuthal angle can be observed. In Fig.^ we show also the role of the evolution 
scale in the Kwieciriski distributions. The QCD evolution embedded in the Kwieciriski evolution 
equations populate larger transverse momenta of partons entering the hard process. This signif- 
icantly increases the initial (nonperturbative) decorrelation in azimuth. For transverse momenta 
of the order of ~ 10 GeV the effect of evolution is of the same order of magnitude as the effect 
due to the nonperturbative physics of hadron confinement. For larger scales of the order of jj} ~ 
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Figure 6: The angular correlations for all four components: gg gg (solid), gg — > qq (dashed) and gq gq 
= qg ^ qg (dash-dotted). The calculation is performed with the Kwiecihski UPDFs and bo= 1 GeV^^ The 
integration is made for jets from the transverse momentum interval: 5 GeV < pi^t,P2,t < 15 GeV and from 
the rapidity interval: -4 <yi,y2 < 4. 

100 GeV^, more adequate for jet production, the initial condition is of minor importance and the 
effect of decorrelation is dominated by the evolution. Asymptotically (infinite scales) there is no 
dependence on the initial condition provided reasonable initial conditions are taken. 
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Figure 7: The azimuthal correlations for the gg — > gg component obtained with the Kwieciiiski UGDFs for 
different values of the nonperturbative parameter b(, and for different evolution scales /x^ = 10 (on line blue), 
100 (on line red) GeV^. The initial distributions (without evolution) are shown for reference by black lines. 

In FigJS] we show the maps for different UGDFs and for gg ggg processes in the broad range 
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Figure 8: Two-dimensional distributions in pit and p2t for KL (left upper), BFKL (right upper), Ivanov- 
Nikolaev (left lower) UGDFs and for the gg — > ggg (right lower). In this calculation -4- < yi,y2 < 4. 

of transverse momenta 5 GeV < pi^t,P2,t < 20 GeV for the RHIC energy W = 200 GeV. In this 
calculation we have not imposed any particular cuts on rapidities. We have not imposed also any 
cut on the transverse momentum of the unobserved third jet in the case of 2 — > 3 calculation. The 
small transverse momenta of the third jet contribute to the sharp ridge along the diagonal pu = P2.t- 
Naturally this is therefore very difficult to distinguish these three-parton states from standard two 
jet events. In principle, the ridge can be eliminated by imposing a cut on the transverse momentum 
of the third (unobserved) parton [Q]. There are also other methods to eliminate the ridge and 
underlying soft processes which is discussed in Ref.[Q]. 

When calculating dijet correlations in the standard NLO (2 —>■ 3) approach we have taken all 
possible dijet combinations. This is different from what is usually taken in experiments where 
correlation between leading jets are studied. In our notation this means p^j < pi^t and p'i^t < P2,t- 
When imposing such extra condition on our NLO calculation we get the dash-dotted curve in Fig.^. 
In this case da/d^^ = for 0_ < 2/3 71. This vanishing of the cross section is of purely kinemat- 
ical origin. Since in the ^, -factorization calculation only two jets are explicit, there is no such an 
effect in this case. This means that the region of ^_ < 2/371 should be useful to test models of 
UGDFs. For completeness in Fig.|I^ we show a two-dimensional plot {pij,P2,t) with imposing 
the leading-jet condition. Suiprisingly the leading-jet condition removes a big part of the two- 
dimensional space. In particular, regions with p2,f > '2-Pi,t and pu > 2p2,f cannot be populated 
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Figure 9: Dijet azimuthal correlations da/d(p^ for the gg gg component and different UGDFs as a 
function of azimuthal angle between the gluonic jets. In this calculation W = 200 GeV and -1 < yi,y2 < 1, 
5 GeV < pit,P2t < 20 GeV. Here the thick-solid line corresponds to the Kwiecinski UGDF, the dashed line 
to the Kharzeev-Levin type of distribution and the dotted line to the BFKL distribution. The two thin solid 
(on line black) lines are for NLO collinear approach without (upper line) and with (lower line) leading jet 
restriction. 




Figure 10: Cross section for the gg ggg component on the {pi.t,P2,t) plane with the condition of leading 
jets (partons). 

via 2 — > 3 subprocess ^ . There are no such limitations for 2 — > 4, 2 ^ 5 and even higher-order 
processes. Therefore measurements in "NLO-forbidden" regions of the {p\^t,P2,t) plane would test 
higher-order terms of the standard collinear pQCD. These are also regions where UGDFs can be 
tested, provided that not too big transverse momenta of jets are taken into the correlation in or- 
der to assure the dominance of gluon-initiated processes. For larger transverse momenta and/or 

'in LO collinear approach the whole plane, except of the diagonal pi , = P2,f' is forbidden. 
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forward/backward rapidities one has to include also quark/antiquark initiated processes via uninte- 
grated quark/antiquark distributions. 

4. Summary 

Motivated by the recent experimental results of hadron-hadron correlations at RHIC we have 
discussed dijet correlations in proton-proton collisions. We have considered and compared results 
obtained with colUnear next-to-leading order approach and leading-order ^^-factorization approach. 

In comparison to recent works in the framework of /^f -factorization approach, we have included 
two new mechanisms based on gq gq and qg qg hard subprocesses. This was done based 
on the Kwiecihski unintegrated parton distributions. We find that the new terms give significant 
contribution at RHIC energies. In general, the results of the fc,-factorization approach depend on 
UGDFs/UPDFs used, i.e. on approximation and assumptions made in their derivation. 

The results obtained in the standard NLO approach depend significantly whether we consider 
correlations of any jets or correlations of only leading jets. In the NLO approach one obtains ^ 
= if < 2/3 71 for leading jets as a result of a kinematical constraint. Similarly ^pf^'^p^^ = if 
Pl,t > 2/?2,r or p2,t > 2pij. 

There is no such a constraint in the -factorization approach which gives a nonvanishing cross 
section at small relative azimuthal angles between leading jets and transverse-momentum asymmet- 
ric configurations. We conclude that in these regions the -factorization approach is a good and 
efficient tool for the description of leading-jet correlations. Rather different results are obtained 
with different UGDFs which opens a possibility to verify them experimentally. Alternatively, the 
NLO-forbidden configurations can be described only by higher-order (NNLO and higher-order) 
terms. We do not need to mention that this is a rather difficult and technically involved computa- 
tion. 

What are consequences for particle-particle correlations measured recently at RHIC requires a 
separate dedicated analysis. Here the so-called leading particles may come both from leading and 
non-leading jets. This requires taking into account the jet fragmentation process. 
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